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I. INTRODUCTION 

Protecting quantum information from damage will be 
essential for future quantum technologies. Thus, devis¬ 
ing a self-correcting quantum memory , capable of stor¬ 
ing qubits for a very long time without active error cor¬ 
rection m , would be a scientific milestone with far- 
reaching implications for the scalability of quantum com¬ 
puting and the security of quantum communication pro¬ 
tocols. 

We define a self-correcting quantum memory as a 
quantum many-body system with a local Hamiltonian 
such that (i) the ground state is degenerate in the limit 
of a large system, (ii) at a sufficiently small nonzero tem¬ 
perature, quantum information stored in the degenerate 
ground state is protected from error for a time which 
grows without bound as the system size increases, (iii) 
both of these features are stable with respect to generic 
small local static perturbations of the Hamiltonian. 

For a quantum memory to be useful, one will also need 
a practical method for storing a quantum state in memory 
and reading it out, but we do not impose that criterion 
in our definition. Rather, we regard (ii) as satisfied if it is 
possible in principle to recover the stored quantum state 
from a thermally fluctuating state. Topologically ordered 
systems have a perturbatively stable ground-state degen¬ 
eracy [3 El, and would meet all three of our criteria if the 
topological order persists at nonzero temperature BUI Hi¬ 
ll- 

Known mathematical models of two-dimensional topo¬ 
logically ordered media (such as the toric code [10]) typ¬ 
ically do not satisfy (ii) mHm- The rate per unit area 
for thermal production of anyons is suppressed by the en¬ 
ergy gap, but does not depend on the system size; once 
produced, thermal anyons can propagate without further 
energy penalty, quickly destroying the quantum informa¬ 
tion stored in the ground state [3 111 IToHTT?] . Three- 
dimensional models have been proposed in which the 
quantum memory time is enhanced because the prop¬ 
agation of excitations is impeded [25H23] . but so far the 


known models of self-correcting quantum memory require 
at least four spatial dimensions mm- (A recently pro¬ 
posed model based on a low-distortion three-dimensional 
embedding of a four-dimensional code [25] has not been 
shown to be perturbatively stable.) For a recent review 
of quantum memories at nonzero temperature, see [27j . 

To suppress errors arising from the propagation of ther¬ 
mally produced anyons, Hamma et al. proposed cou¬ 
pling the toric code to massless bosonic fields, inducing 
long-range attractive interactions among the anyons |28j ; 
however, their scheme requires the coupling strength to 
diverge with the system size. Pedrocchi et al. pro¬ 
posed an alternative scheme, in which coupling anyons 
to three-dimensional massless bosons induces an infrared 
divergent renormalization of the anyon chemical poten¬ 
tial, forcing the rate for thermal anyon production to 
approach zero as the system size grows [2M2T] . Another 
suggestion is that disorder in a two-dimensional topolog¬ 
ical medium could localize the anyons and hence extend 
the memory time [321 i .'TTSl . 

In this paper, we emphasize a general issue confronting 
schemes for stabilizing a two-dimensional quantum mem¬ 
ory by coupling to massless bosons — the criterion (iii) 
typically fails because generic small local perturbations 
introduce a mass gap for the bosons and hence compro¬ 
mise the scheme. We consider the model of Pedrocchi 
et al. [31] as an instructive example, but our arguments 
apply more generally, and we also assess some alternative 
proposals for enhancing memory times. We offer no rig¬ 
orous arguments excluding self-correcting memory in two 
or three dimensions, but we hope our discussion helps to 
clarify the question. Nor do we mean to deny the possi¬ 
bility that clever engineering of interactions or disorder 
in anyonic systems could substantially enhance the mem¬ 
ories times of anyonic systems; our focus is on whether 
arbitrarily long memory times are achievable as a matter 
of principle. 

This paper is organized as follows. In Section [n] we 
discuss how Hamiltonian perturbations and thermal fluc¬ 
tuations can corrupt the information stored in a quantum 
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memory, and in Section III we recall some proposals for 
enhancing the storage time. In Section IV we discuss 
the idea of stabilizing a two-dimensional topologically 
ordered system at nonzero temperature by introducing 
long-range interactions between anyons, using the “toric- 
boson model” introduced in [3lj as an explicit example, 
and in particular we explain how to recover the conclu¬ 
sions of [3TJ using an effective field theory description of 
the system. In Section [V] and [VlJ we argue that the in¬ 
finite energy barrier of the toric-boson model, and hence 
its thermal stability, are unstable with respect to generic 
perturbations of the system’s microscopic Hamiltonian. 
Section VII contains some concluding remarks and open 
questions regarding self-correcting quantum memory. 


II. LOGICAL ERRORS IN A QUANTUM 
MEMORY 

To clarify the concept of self-correcting quantum mem¬ 
ory, we recall in this Section some ways for a quantum 
memory to fail, either due to imperfect control of the 
Hamiltonian or due to thermal fluctuations. 


A. Spectral instability 

Desideratum (iii) of our definition of a self-correcting 
memory requires the degenerate ground space to be 
stable with respect to small local perturbations of the 
Hamiltonian. If adding a small local perturbation can 
lift the degeneracy, then the environment might couple to 
this local operator, producing fluctuations in the relative 
phases of energy eigenstates, hence driving decoherence 
in the energy-eigenstate basis. This failure mechanism 
relies only on the properties of the nearly-degenerate 
ground space, not the properties of the excitations sup¬ 
ported by the system. 

Alternatively, the memory could fail because a weak lo¬ 
cal perturbation destroys the energy gap separating the 
ground space from the rest of the energy-eigenstate spec¬ 
trum. When the energy gap collapses, the structure of 
the ground space may undergo a qualitative change, al¬ 
lowing the information stored in the ground state to be 
corrupted. 

We use the term “spectral instability” to refer to either 
of these failure modes. Sufficient conditions for spectral 
stability, based on topological order, were identified in |6] 
and then generalized in |34j . 


B. Thermal relaxation 

Desideratum (ii) of our definition requires the stored 
quantum information to be well protected from error at a 
sufficiently small nonzero temperature. Ideally, we would 
assess whether (ii) is satisfied by modeling the thermal 
environment, simulating how the stored information is 


affected by thermal fluctuations, and computing the fi¬ 
delity of a suitably chosen decoding map. Since this pro¬ 
gram is quite difficult to carry out, and in any case the 
results depend on the details of the environment and de¬ 
coding procedure, we usually assess (ii) using simpler cri¬ 
teria. 


1. Energy and free energy barriers 

One way to simplify the story is to imagine that the 
environment applies a sequence of local quantum opera¬ 
tions to the quantum memory. After many such steps, 
we envision cooling the memory, so that it relaxes back 
to the ground-state space, and then we compare this final 
state to the initial state deposited in the memory. We say 
the error sequence is “harmful” if the final state deviates 
significantly from the initial state. 

Each error sequence has an energy height, the max¬ 
imal value of the expectation value of the Hamiltonian 
attained during the sequence, and we say that the mem¬ 
ory’s energy barrier is the lowest energy height for any 
harmful error sequence. Naively, a high energy barrier 
means a long memory time, because only rare thermal 
fluctuations will be able to surmount the barrier and 
damage the stored quantum state. 

An Ising ferromagnetic in two or more spatial dimen¬ 
sions is a thermally stable classical memory storing one 
bit of information in the sign of its magnetization; for the 
magnetization to flip, thermal fluctuations must produce 
a large droplet of flipped spins with size comparable to 
the system size, surmounting an energy barrier scaling 
like the droplet’s surface area. The divergent memory 
time, though, is not perturbatively stable. A small mag¬ 
netic field favors one sign of the magnetization; hence 
bubbles of flipped spins, arising as thermal fluctuations 
in the disfavored phase at a constant rate per unit time 
and volume, destroy the stored classical information in a 
constant time. The three-dimensional toric code, on the 
other hand, is an example of a self-correcting classical 
memory, as the one bit detected by its membrane-like 
logical operator is robust against small local perturba¬ 
tions of the Hamiltonian. 

As in the classical case, to ensure a long storage time 
for a quantum memory we would like the energy barrier 
to diverge as the system size increases. However, a diver¬ 
gent energy barrier is not by itself sufficient for thermal 
stability. Haah’s three-dimensional cubic code [20] has a 
logarithmically divergent barrier and Michnicki’s welded 
code [25, 53] has a polynomially divergent barrier, yet 
in both cases the memory time is bounded above by a 
constant independent of the system size at any nonzero 
temperature EH [27]. 

These examples caution us that even if harmful error 
sequences are strongly Boltzmann suppressed, the sup¬ 
pression can be overwhelmed if the harmful sequences 
have high entropy. We should consider not just the en¬ 
ergy cost, but rather the free energy cost of encoded er- 
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rors. In the case of the cubic code, for example, the 
Boltzmann suppression dominates when the system size 
is sufficiently small; hence the logarithmically growing 
energy barrier implies that the memory time r increases 
with system size L as r ~ exp (c/3 log L ), where c is a con¬ 
stant and (3 is the inverse temperature. However once L 
is larger than a critical system size L* ~ e c entropy be¬ 
comes dominant, and the memory time actually declines 
with increasing system size. The memory time is opti¬ 


mized for L ~ L *, yielding r < e c @ for some constant c; 
criterion (ii) for self correction is not satisfied. 

Conceivably, entropic effects might enhance the mem¬ 
ory time under some circumstances, perhaps allowing the 
memory to be self correcting even though the energy bar¬ 
rier is finite, as suggested in [IS]. A model in which en¬ 
tropy suppresses ha rmful error sequences was proposed 
in [55] (see also Sec. Ill B l. 


2. Anyon propagation 


R requires energy scaling like logi?; thus thermal par¬ 
ticle propagation is suppressed and the memory has a 
logarithmically growing energy barrier. Nevertheless, as 
noted in Sec. mm the free energy cost of a logical er¬ 
ror becomes independent of system size for a sufficiently 
large system at a fixed temperature. 

In the three-dimensional welded code [22, (23], v > 2 
three-dimensional toric code subblocks are sewn together 
at each weld. Anyons propagate freely within each sub¬ 
block, but for an anyon to pass through a weld from 
one subblock to another, at least v — 2 additional anyons 
must be created. Therefore the energy barrier is high 
if the code has many such subblocks. However, logical 
errors can arise at nonzero temperature from a process 
in which v thermal anyons propagating in distinct sub¬ 
blocks meet at a weld and fuse to the vacuum [27J. This 
process is entropically enhanced because the fusion can 
occur at any point along the weld; the resulting logical 
error rate is independent of system size and therefore the 
code is not self correcting. 


Typical topologically ordered systems in two dimen¬ 
sions support anyons — pointlike quasiparticle excita¬ 
tions obeying exotic statistics. If the system resides on a 
torus, a logical error occurs when a pair of anyons is cre¬ 
ated, and then one member of the pair propagates around 
a homologically nontrivial cycle of the torus before rean¬ 
nihilating with its partner. Thus the energy barrier is 
the (constant) energy cost of creating and separating a 
pair of anyons. 

At nonzero temperature, thermal fluctuations produce 
a nonvanishing density of anyons per unit area. We can 
imagine decoding the state by bringing pairs of anyons 
together to annihilate, returning the system to a ground 
state. Starting with an initial noisy encoded state, we 
may allow the system to evolve for a while in contact 
with a heat bath, before applying the decoding map to 
the final configuration. By combining the anyon world 
lines arising from the decoding of the initial state, the 
decoding of the final state, and the thermally induced 
anyon propagation, creation, and annihilation between 
the initial and final time, we obtain a cycle, a set of 
closed loops. If this cycle is homologically nontrivial, a 
logical (encoded) error occurs. 

When the anyons propagate a distance comparable to 
the typical (constant) distance between particles in the 
thermal anyon gas, a logical error will occur with nonneg- 
ligible probability. Therefore, we do not normally expect 
a two-dimensional quantum memory to be self correct¬ 
ing, unless supplemented by an additional mechanism to 
impede the propagation or creation of anyons. 

The three-dimensional version of the toric code also has 
pointlike “anyon” excitations, and a logical error occurs 
when an anyon propagates across the system. Therefore 
this code, too, has a constant energy barrier and is not 
self correcting. The cubic code [20] has pointlike exci¬ 
tations, but reaching a configuration in which one par¬ 
ticle is separated from all others by a distance at least 


III. IMPROVING THE STORAGE TIME 

Having discussed how a quantum memory can fail, we 
now briefly consider some possible ways to improve the 
stability of a quantum memory and hence enhance its 
storage time. 


A. Anyon localization 

Even for a topologically ordered system with a spectral 
gap, in which small local perturbations split the ground 
state degeneracy by an amount which is exponentially 
suppressed for a large system size, a quantum memory 
may fail due to small local perturbations of the Hamilto¬ 
nian. The problem is that, although information encoded 
in the ground space may be well protected in principle, 
if the perturbation is unknown it may be difficult to pre¬ 
pare precisely the ground state of the perturbed system. 

For example, if the ideal system has the toric code 
Hamiltonian, we might write to the memory by prepar¬ 
ing one of the toric code ground states. But if the Hamil¬ 
tonian is slightly deformed by a weak perturbation, the 
ground state of the ideal Hamiltonian is a coherent su¬ 
perposition of excited states of the perturbed Hamilto¬ 
nian; this superposition dephases quickly, destroying the 
encoded information in a constant time independent of 
system size 1.3(11 . If the perturbation were exactly known, 
then in principle we could decode the state by running the 
Hamiltonian evolution backward in time. But if (more 
realistically) the perturbation is unknown, then an un- 
correctable logical error is likely to occur. This error 
mechanism operates even at zero temperature; the er¬ 
ror arises not from thermal fluctuations but rather from 
imperfect control of the system’s Hamiltonian. 
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One way to evade this problem is to respond to the de¬ 
formation of the Hamiltonian by adjusting the procedure 
for encoding and decoding the state. Under a generic per¬ 
turbation, the static pointlike anyons of the exact toric 
code become propagating dressed anyons with nonzero 
size. Using a physical device that detects these dressed 
anyons, we can prepare a ground state of the perturbed 
system, which will store quantum information reliably at 
zero temperature. In principle we could read out the en¬ 
coded state at a later time using dressed string operators. 

Alternatively, these zero-temperature errors can be 
suppressed by introducing disorder in the system Hamil¬ 
tonian, ensuring that in low-lying excited states the 
anyons are spatially localized rather than freely prop¬ 
agating [33]. In fact, localization may suffice to make 
the memory time exponentially large in the system size, 
when the errors are due to unknown Hamiltonian pertur¬ 
bations [SUETj. However, at any nonzero temperature, 
single-particle localization does not prevent anyons from 
diffusing at a nonzero rate due to thermally activated 
hopping through the disordered system. Thus at nonzero 
temperature, though localization may enhance the mem¬ 
ory time by a constant factor, we still expect the logical 
errors to occur after a constant time independent of sys¬ 
tem size — the memory is not self correcting. 

We note that a disordered topologically ordered sys¬ 
tem may exhibit many-body localization, in which anyons 
are spatially localized even in highly excited states of 
the disordered Hamiltonian [33] [35]. However, many- 
body localization is a property of a perfectly isolated 
Hamiltonian system. Even a many-body-localized two- 
dimensional topological phase would fail to be self cor¬ 
recting, because anyons would still be able to diffuse 
when the system is in contact with a thermal bath. 

For the case of a two-dimensional commuting Hamilto¬ 
nian like the disordered toric code model, the finding that 
topological order is destroyed at nonzero temperature can 
be formulated rigorously. Hastings [7] showed that the 
Gibbs state of a (not necessarily translation-invariant) 
2D commuting local Hamiltonian can be (approximately) 
transformed to a product state by a constant-depth local 
quantum circuit. To analyze the memory time, we can 
draw guidance from studies of the disordered ID Ising 
model, whose stability properties are similar to those of 
the 2D toric code. For a one-dimensional classical spin 
system with Hamiltonian 

Hlsing = ^ ^ (-1 ( 1 ) 

i 

where the nearest-neighbor couplings {Ji} are positive 
but otherwise arbitrary, the memory time r is indepen¬ 
dent of system size, scaling like r ~ exp(2/3J) where 
J = ( J m in + Jmax)/ 2 is the average of the minimal and 
maximal values of the coupling [40] . Thus the mem¬ 
ory time is not significantly enhanced in the regime 
Jmax Jmin where defects are strongly localized. The 
conclusion r = 0(1) also applies to the (disordered) toric 


code, subject to a master equation whose jump operators 
are Pauli operators; in that case the jumps transform en¬ 
ergy eigenstates to energy eigenstates, so that the con¬ 
vergence to a fixed point can be analyzed using classical 
methods mi- 

B. Entropically suppressed propagation 

Brown et al. described an anyon model in which anyon 
propagation is impeded by entropic considerations [35]. 
In their model, a two-dimensional topologically ordered 
medium is divided into domains separated by defect lines. 
Each domain supports both heavy and light anyons, 
where the decay of a heavy anyon into two light anyons 
is allowed by both kinematics and the fusion rules. Fur¬ 
thermore, upon crossing a defect line, a heavy anyons is 
transformed into a light one, and a light one is trans¬ 
formed into a heavy one. 

Propagation of light anyons across defect lines is sup¬ 
pressed at low temperature, because an energy boost is 
needed to cross to the other side. Furthermore, it is en¬ 
tropically favored for heavy anyons to decay to pairs of 
light ones. Therefore, anyons cannot easily travel across 
the system and the memory time is correspondingly en¬ 
hanced. Numerical evidence suggests that if the system 
size is properly optimized, the memory time r scales with 
the inverse temperature fj as r ~ e c/3 , just as in the 
three-dimensional cubic code, even though the energy 
barrier is constant rather than logarithmically growing. 
This is a notable improvement compared to the scaling 
r ~ e c ^ in the toric code. Though not self correcting, this 
system illustrates how in some cases favorable entropic 
effects can significantly enhance the memory time. 

C. Long-range anyon interactions 

Since logical errors in a two-dimensional topologically 
ordered quantum memory arise from anyon propagation, 
any mechanism that impedes anyon propagation can in 
principle enhance the memory time. Aside from localiza¬ 
tion due to disorder as discussed above, there is another 
mechanism that could limit how far anyons can travel 
— attractive interactions between anyons. The interac¬ 
tions may impose a hefty energy penalty on harmful er¬ 
ror sequences, encouraging confinement of the anyons to 
small dilute clusters which do not produce logical errors. 
For the system to be truly self correcting, though, these 
interactions must have infinite range. Otherwise only a 
constant amount of energy would be needed for an anyon 
to break away from a cluster, leaving it free to propagate 
unmolested by attractive forces applied by other anyons, 
thus allowing logical errors to occur at a constant rate at 
any nonzero temperature. 

For a system with a local Hamiltonian, long-range 
forces arise from the exchange of gapless excitations. 
Hamma et al. [28] proposed a model in which the anyon 
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density provides a source for a massless scalar field. 
Scalar exchange produces long-range attractive interac¬ 
tions between anyons, but each anyon also has a logarith¬ 
mically infrared divergent negative self energy, rendering 
the system unstable. The divergent self-energy can be 
cancelled by introducing other fields with carefully tuned 
couplings, but then the Hamiltonian contains terms with 
infrared divergent coefficients. 

Improving on j2<§| . Pedrocchi et al. proposed an al¬ 
ternative model with a local Hamiltonian and finite cou¬ 
plings, in which the anyonic chemical potential is infrared 
divergent [31]. Thus the energy barrier, the cost of 
creating a distantly separated anyon pair, diverges with 
system size, as does the memory time. We review this 
scheme in Sec. IV It works, if the Hamiltonian is care¬ 
fully tuned. But as we explain in Sec. [V] it does not 
meet our criteria for self correction, because the diver¬ 
gent energy barrier is unstable with respect to generic 
small perturbations of the local Hamiltonian. 

We also note that schemes for “active” rather than 
“passive” error correction have been proposed, inspired 
by the idea that long-range attractive interactions might 
impede the propagation of diffusing anyons [31]. In a 
scheme recently discussed in [32] . the anyon motion is 
guided by a classical field, which is regularly updated us¬ 
ing a local evolution rule as the anyons propagate. Such 
active error correction schemes, though effective and in¬ 
teresting, do not meet our criteria for self correction for 
several reasons. For one thing, to achieve an arbitrar¬ 
ily low logical error rate, the scheme requires a “speed 
of light” (rate of information propagation in the classi¬ 
cal field) which diverges relative to the “speed of sound” 
(rate of anyon motion), a requirement not satisfied by a 
system with a local Hamiltonian and bounded couplings. 
Furthermore, the dissipative dynamics of the field which 
controls the anyons may not be realizable by any Hamil¬ 
tonian system in contact with a thermal bath. 


IV. TORIC CODE COUPLED TO A SCALAR 
FIELD 


A. Toric code 


To prepare for our discussion of models with long- 
range anyon interactions, let’s first recall Wen’s formula¬ 
tion [35] of Kitaev’s toric code model m- Consider an 
L x L square lattice with one qubit per site, where L is 
an even integer. The Hamiltonian is 

H-Toric = J ^ ' If /-■ W r = i Z r ^2-A r ^Z r 4 j ( 2 ) 
r 

where r labels lattice plaquettes and W r acts on the 
four qubits contained in plaquette r (X and Z denote 
Pauli matrices). On the torus ( i.e. a square with peri¬ 
odic boundary conditions), the model has four degenerate 
ground states which can encode two qubits. 


The energy cost of creating an anyon pair is 4J, and 
the anyons are noninteracting. A logical error occurs if a 
pair is created and one anyon propagates around a non¬ 
trivial cycle of the torus before the anyons reannhihilate. 
Therefore the (constant) energy barrier is 4J. The typ¬ 
ical separation between particles in the thermal anyon 
gas, and hence also the memory time, scale like the in¬ 
verse Boltzmann factor e 2/3J . 


B. Coupling the toric code to massless bosons 


Now, following Pedrocchi et al. [31], we consider cou¬ 
pling the two-dimensional toric code to an auxiliary 
three-dimensional bosonic system which induces long- 
range attractive forces among the plaquette variables. 
The Hamiltonian proposed by [31! is 


H — A ^ [ W r (a r T q£) T Hbosoni 


Hboson — ^ ^ ^ 


qOqjO'q, 


(3) 


a * = ^E a * e 


iqs 


Here the index r labels the sites of the two-dimensional 
topologically ordered medium and s labels the N sites 
of a three-dimensional lattice which contains the two- 
dimensional lattice; a r and a\ are the annihilation and 
creation operators for bosonic modes residing at the loca¬ 
tion of plaquette operator W r , and q labels the discrete 
Fourier modes for a bosonic field defined on the three- 
dimensional lattice. We assume e q > 0, because a bosonic 
mode with negative e q would be unstable. The precise 
behavior of e q for large q does not matter for analyzing 
the long-distance physics of the model, but the behav¬ 
ior of e q for small q does matter. We assume a gapless 
spectrum for the Fourier modes with e Q ~ q 2 for small q. 
We will call the Hamiltonian in Eq. (pi the “toric-boson 
model.” 

Coupling anyons to the three-dimensional massless bo¬ 
son field produces an attractive force between anyons, 
with potential energy inversely proportional to separa¬ 
tion. After a polaron transformation (i.e., completing 
the square) the Hamiltonian may be rewritten as [31] 


H = -J2 JryW r W r , + ^ e q a\a q , (4) 


where 


A 2 ^ 1 


Jry - N E , 


,iq(r-r') 


A 2 


r — r 


(5) 


here N is the number of lattice sites in the three- 
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dimensional bosonic system, 

5, = a, --^= (6) 

and the conclusion J r ^ r j ~ A 2 /\r — r'\ applies for asymp¬ 
totically large separation. Note that d q and cA q both obey 
ordinary bosonic commutation relations because the W r ’ s 
residing at distinct plaquettes are commuting operators. 
The long-range behavior of the potential energy, scaling 
like ~ 1/7? where R is the separation, arises from the pole 
in 1 /c q , which dominates the Fourier integral for large R. 

Because the Hamiltonian in Eq. Q is the sum of an 
anyonic part Ha and bosonic part Hb, which are com¬ 
pletely decoupled, the energy eigenstates can be written 
as I ip) = I ipqubit) <8> \(t>b oson ) •> where \ip qu bit) is an eigen¬ 
state of Ha and 1 4>boson) is an eigenstate of Hb- Since 
J r y > 0, the interactions between plaquette variables 
in Eq. Q are “ferromagnetic,” favoring alignment. Ha 
has a Z 2 global symmetry, under which W r changes sign, 
which is spontaneously broken at zero temperature - 
there is one ground state sector with W r = 1 for all r 
and another sector with W r = — 1 for all r, where both 
sectors have a topological degeneracy as well. The spon¬ 
taneous breakdown of the global Z 2 symmetry persists 
at sufficiently small nonzero temperature, and we may 
assume without loss of generality that the thermal ex¬ 
pectation value (W r ) is positive . 

As emphasized in [3T], the crucial property of Eq. © 
is that the energy cost of creating a single isolated anyon 
is infrared divergent, i.e., increases with the system size. 
To compute this energy cost, we follow .31] and suppose 
that the bosonic bath is at any moment in a Gibbs state 
with respect to the effective Hamiltonian (a\H\a), where 
|a) is the instantaneous anyonic state, i.e., an eigenstate 
of all operators W r . In that case, this energy cost of cre¬ 
ating a single isolated anyon only depends on the anyonic 
part Ha of the Hamiltonian in Eq. Q (see m for the 
detailed calculation). The divergence arises because each 
plaquette interacts with many other distantly separated 
plaquettes on the two-dimensional lattice; flipping the 
sign of one W r relative to all the others increases the 
potential energy by an amount 

~0(L). (7) 

r' 

This effective chemical potential /j, is O(L) because we 
sum the 1/R interaction energy over the Lx L lattice. In 
contrast to the standard toric code, where the energy of 
an isolated anyon is 2 J and the thermal density of anyons 
at inverse temperature /3 scales like e 2/3J , in the model 
Eq. 0 the thermal density scales like e c/3i (where c is a 
constant), and correspondingly the memory time grows 
exponentially with L at any nonzero temperature. The 
infrared divergent chemical potential and the resulting 
self-correcting behavior of Eq. © were the main con¬ 
clusions of |5T] . In fact, the infinite chemical potential 


means that the spontaneous breakdown of the Z 2 global 
symmetry, and hence the self-correcting behavior, hold 
at arbitrarily high temperature. 


C. Derivation using effective field theory 

Because the conclusion y. = O(L) follows from only 
the long-distance properties of the Hamiltonian Eq. Q, 
it can be recovered using an effective field theory descrip¬ 
tion which captures the infrared behavior of the model. 
The effective Hamiltonian becomes 

Heft = \ J d 3 x (V0) 2 - A J d 2 yw{y)(j){y) H-, (8) 

where cj)(x ) is a suitably normalized scalar field. The 
| (V</>) term in the Hamiltonian density arises from the 
dispersion relation e q ~ q 2 for small g; higher-order terms 
in q correspond to terms with more derivatives in Heft , 
which are not shown because they produce only small 
corrections for long-wavelength modes of the scalar field. 
The source w(y) describes the anyon configuration. To be 
compatible with the conventions of Sec. |IVB[ we assume 
w(y) takes the positive value Wq > 0 in the ground state 
of the two-dimensional topological medium, and has the 
value w(y) = wo — 2n(y) when anyons are present, where 
n(y) is a suitably smoothed and normalized anyon den¬ 
sity. We assume w(y) =0 outside the two-dimensional 
medium. 

The energy is minimized by a solution to the static 
classical field equation, which is Poisson’s equation 

\I 2 (j){x) = —Aw{ x). (9) 

The solution is 

4>(x) = A f d 2 y ( 3 ^ 2 ) w(y) + (j) 0 (x) (10) 

where <f>o(x) is a solution to the homogeneous field equa¬ 
tion, and the Green function • ^. 2 is determined by im¬ 
posing appropriate boundary conditions. We assume 
that the two-dimensional anyonic source for the three- 
dimensional scalar field </> has bounded support, and that 
(j) vanishes at spatial infinity. 

The boundary condition justifies an integration by 
parts with no surface term, and using the field equation 
we may obtain an alternative expression for the minimal 
energy of the field configuration in the presence of an 
anyon source: 

E[w\ = - j d 3 x\^ 2 ct>-A J d 2 yw{y)c/)(y) 

= | J d 3 xw(x)cj){x) - A J d 2 yw{y)<t>{y) 

= (ID 
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with the integral now supported only on the two- 
dimensional surface where the source field w is nonzero. 


Inserting Eq. (10) into this expression, one finds 


E[w\ = -4/ [ d 2 xd 2 yw(x) ( 3 ^ 2 ) w(y). (12) 

J ' ' x,y 


The boundary condition determines the Green function 
to be 



1 

47t \x — y\' 


(13) 


and so we obtain 


E[w\ = -\ 


d 2 xd 2 yw(x)J(x,y)w(y), (14) 


where 


J{x ’ v) = <15) 

as in Eq. ©• 

To compute the effective anyon chemical potential, we 
write w(y) = wq — 2 n(y), finding that E[w] contains a 
term linear in the anyon density 

I*** w " w (16) 

where 

d’i.x) = 2wq J d 2 y J(x, y). (17) 


Assuming for simplicity that x is the center of a disk with 
radius A/2, we obtain 

/»27t /»L/2 

y = 2 w 0 A 2 / dO rdrH- = L , (18) 

Jo Jo 

recovering the linear divergence found in Eq. (J7|. 


V. PERTURBATIVE INSTABILITY OF 
MASSLESS BOSONS 


Having seen that the model in Eq. ([ 3 ]) has an infrared 
divergent energy barrier, and hence an exponentially 
increasing quantum storage time at sufficiently small 
nonzero temperature, we next consider the stability of 
the energy barrier with respect to small perturbations of 
the system’s local Hamiltonian. For assessing perturba¬ 
tive stability, it is important to note that the divergent 
anyon chemical potential arises because the response of 
the scalar field to the anyon vacuum is infrared divergent. 
A small perturbation of the Hamiltonian can tame this 
divergence, rendering the energy barrier finite. 


The bosonic ground state \4>b OS on) of Eq. Q satisfies 

dq\4>boson ) = 0 ( 19 ) 

for all q , which implies 

K>0 ~ ^ E r E eMr '~ r) = °( L ) (2°) 

q H r' 

where (a r )0 denotes the ground-state expectation value, 
and we have substituted W r > = 1. The bosonic mode 
a r has an unbounded occupation number in the anyon 
ground state, which is why the cost of creating an isolated 
anyon at site r is divergent. 

Additional terms in the bosonic Hamiltonian, which 
are generically present, will prevent the occupation num¬ 
ber from diverging. Suppose, for example, the Hamilto¬ 
nian contains a term 

V=e E ,lr=£ E fl 9 S 9> ( 21 ) 

r q 


where n r = a' r a r is the number operator for the mode lo¬ 
calized at site r. We may assume e > 0 because otherwise 
the bosonic field would be unstable, an indication that we 
are expanding around the wrong bosonic vacuum. Now 
the perturbation V, which disfavors a large occupation 
number, competes with the Hamiltonian Eq. Q, and the 
ground-state expectation value of a r takes the finite value 


{Q’r} 0 
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1 

e + Eg 


^y«( r '- r ) = 0(e -1 / 2 ), 

r' 


( 22 ) 


because the perturbation removes the pole at q = 0. 

For this particular perturbation the bosonic Hamilto¬ 
nian is quadratic and therefore exactly solvable, but the 
conclusion that a generic perturbation removes the in¬ 
frared divergence holds more generally. In the effective 
field theory language, the solution to the static field equa¬ 
tion Eq. (10) becomes 


4>{x) = Afd 2 y ( 3 ^ 2 ) w(y) 

\ / x,y 

= Aw 0 fd 2 yi ^^=O(L) (23) 


in the anyonic vacuum such that w(y ) = Wq in the two- 
dimensional topological medium and w(y) = 0 elsewhere. 
The infrared divergent behavior of 4>(x) arises because 
the scalar field is exactly massless in Eq. ([8|. Gener¬ 
ically a nonzero mass term will be present, perturbing 
the Hamiltonian according to 


Heft -a Heft + 



\m 2 (j) 2 . 


(24) 


We may assume m 2 > 0, because otherwise the scalar 
field would be unstable, indicating that we are expanding 
around the wrong bosonic vacuum. The static vacuum 









solution becomes 


(j){x) = Awo J d 2 y ( _ v i +m 2 ^ xy = 0(m x ); (25) 

now the y integral converges because the Green function 
decays like e~ m ' x ~ v ' for large separation. 

The effective held theory approach is useful for ana¬ 
lyzing the energy barrier of a quantum memory because 
whether the energy barrier diverges or not hinges on the 
behavior of the model in the far infrared, where the held 
theory approximation should be reliable. The degrees 
of freedom described by the held theory may be emer¬ 
gent, not corresponding in any very direct with the mi¬ 
croscopic degrees of freedom of the underlying Hamil¬ 
tonian model. But just as the divergent energy barrier 
in the model Eq. ([ 8 ]) is realized only when the free pa¬ 
rameter mr of the effective held theory is tuned to zero, 
we may anticipate that a divergent energy barrier arises 
from the exchange of massless bosons only on (at best) a 
codimension-one surface in the space of all Hamiltonian 
models. A generic perturbation moves the Hamiltonian 
off this surface, destroying the inhnite energy barrier. To 
evade this conclusion, we need models in which bosons 
with the appropriate couplings remain gapless when the 
microscopic Hamiltonian is deformed. 

Not just the inhnite energy barrier but also the topo¬ 
logical order of the two-dimensional medium might be 
threatened by small perturbations of the Hamiltonian. 
In addition to coupling to plaquette variables, the scalar 
held could couple to individual qubits; for example the 
perturbed Hamiltonian might contain terms such as 

V = e r 'y ' (a r + a\.)(X r i + X r 2 + X r 3 + X r ^) t (26) 


If in the ground state the expectation value of the held 


scales like (a r ) 0 ~ m 1 ~ e 1 / 2 as in Eq. (22) and (25), 
then an effective magnetic held 


B 


(27) 


is applied to the qubits. Although enhanced by the large 
strength of the held, this magnetic held is still pertur- 
batively small compared to the energy cost ~ 1/y/e of 
creating an isolated anyon, so it is plausible that the 
topological order would survive. In any case, we expect 
the topological order to persist within a region of hnite 
volume in the space of Hamiltonians. 


ceptions to this rule, and we should consider whether 
these exceptions can be exploited by quantum engineers. 
In particular, the masslessness of Goldstone bosons and 
gauge bosons follows from general symmetry principles. 


A. Goldstone bosons 


A spontaneously broken exact continuous symmetry is 
always associated with an exactly massless Goldstone bo¬ 
son, whose masslessness is preserved by any perturbation 
that respects the symmetry (and also maintains its spon¬ 
taneous breakdown). If <j>(: r) is a Goldstone held, then 
the symmetry acts according to 


(j){x) 1 ^ </>( x) + constant, 


(28) 


disallowing a mass term ^m 2 <p 2 in the effective Hamilto¬ 
nian. 

The symmetry not only forbids a mass term for the 
Goldstone boson, but also constrains its coupling to 
other helds. For example, a term —Acf>(x)w{x) describ¬ 
ing the coupling of the Goldstone held to the anyonic 
density w(x) is not allowed, but a derivative coupling 
—AV(f>(x)w(x) is compatible with the symmetry, and can 
be reexpressed as A(j>{x)S7w(x) after integrating by parts. 
(This term violates rotational invariance, but could be 
allowed in an anisotropic system.) Then, arguing as in 
Sec. |IV C we hnd in place of Eq. (12 1 that the energy 
of the bosonic ground state in the presence of a source 
becomes (in D spatial dimensions) 

E[w\ = J d D xd D yXw(x) (^ 2 ) 'Vw(y) 

\ / x,y 

= f d D xw(x) 2 , (29) 


the infrared-hnite integral of a local source-dependent 
energy density. Because they are derivatively coupled, 
the exchange of Goldstone bosons generates short-range 
contact interactions between anyons rather than a long- 
range forces. 

Perturbations that explicitly break the continuous 
symmetry would gap out the Goldstone bosons in any 
case. But even if we are willing to impose a symmetry 
that protects the massless Goldstone bosons, there would 
be no resulting divergent energy barrier to suppress log¬ 
ical errors. 


B. Gauge bosons 


VI. STABILIZING MASSLESS BOSONS 

The scheme of Pedrocchi et al. [3TJ for stabilizing a 
quantum memory requires massless bosonic helds, and 
we have argued that under generic conditions bosonic 
helds acquire mass unless the Hamiltonian is tuned to 
a codimension-one surface. But there are important ex- 


In some cases, exactly massless deconhned t/(l) gauge 
bosons can emerge from an underlying spin model with 
no fundamental gauge symmetry. From an effective held 
theory viewpoint, since the gauge symmetry is not an 
exact feature of the underlying microscopic Hamiltonian, 
nothing seems to forbid adding a gauge non-invariant 
term to the Hamiltonian of the long-distance effective 
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field theory, which would gap out the gauge bosons. But 
if the microscopic local Hamiltonian energetically favors 
a sector of the theory which satisfies gauge constraints, 
these constraints could be robust against generic small 
perturbations of the microscopic Hamiltonian [U- A 
model with such features can support a stable “Coulomb 
phase,” where a “photon” remains exactly massless with¬ 
out any need to impose symmetries or carefully tune the 
Hamiltonian. 

Emergent photons can mediate long-range forces be¬ 
tween charged particles, which are attractive for parti¬ 
cles of opposite charge, and repulsive for particles of like 
charge. Chesi et al. suggested that long-range repulsive 
interactions among anyons in a two-dimensional medium 
could produce a divergent free energy barrier, because for 
a nonzero anyon density the cost of introducing an addi¬ 
tional isolated anyon would diverge with system size j29| ; 
thus the anyon density should approach zero in the ther¬ 
modynamic limit. 

However, anyon models in which gauge-mediated in¬ 
teractions are repulsive for each anyon pair do not seem 
to be self consistent. Because the photon couples to a 
locally conserved charge, anyon pairs which can be lo¬ 
cally created from the vacuum must have opposite “elec¬ 
tric charges.” At nonzero temperature, anyons coupled 
to emergent photons would form a plasma with charge 
screening and no divergent energy barrier 

It may be possible in principle for anyons to carry 
dipole moments which couple to emergent photons. In 
that case, if the gauge field is three-dimensional there 
would be a long-range interaction with potential energy 
falling of with distance R like 1/R 3 . Even if in a suitably 
engineered Hamiltonian the 3D gauge field couples to the 
plaquette variables of the 2D toric code, these interac¬ 
tions decay too rapidly to produce an infrared divergent 
energy barrier. 

For the 3D toric code, a logarithmically divergent 
chemical potential /r ~ logL due to dipolar interactions 
is conceivable. It’s not clear, however, how to arrange the 
needed coupling between the topological medium and the 
gauge theory; see Sec. VII E| for further discussion. 


stability of the long-range interactions, we may consider 
the consequences of long-range forces between “anyons” 
in D dimensions described by Eq. 0 , where 


(30) 


For 0 < a < D, the cost of creating an isolated anyonic 
excitation diverges with the system size, so the model has 
a quantum memory time which increases rapidly with 
system size, even at arbitrarily high temperature. The 
precise scaling of the lifetime of the memory with system 
size is discussed in [T5]. For a > D, energy cost remains 
finite, and the memory time is bounded above by a con¬ 
stant at any nonzero temperature. The case a = D is 
marginal, with a logarithmically divergent energy bar¬ 
rier. 

By similar reasoning, the one-dimensional Ising model 
with 0 < a < 1 is a good classical memory at arbitrarily 
high temperature. It is also known that this system has a 
finite-temperature phase transition for 1 < a < 2, and is 
magnetically disordered at any nonzero temperature for 

a> 2®0Z|. 


VII. DISCUSSION 


The main thrust of this paper is contained in Sec. m 
VI where we have described schemes for stabilizing a 
two-dimensional topologically ordered quantum memory 
at nonzero temperature using long-range interactions, 
and have argued that such schemes have a significant and 
generic drawback — the Hamiltonian needs to be pre¬ 
cisely tuned for the memory time to scale favorably with 
the system size. Our arguments are incomplete, and leave 
many open questions, including the important and en¬ 
ticing question whether self-correcting quantum memory 
is achievable in a suitably engineered three-dimensional 
system. This concluding section contains a few questions 
and remarks which may help to guide future research on 
self-correcting quantum memory. 


C. Other long-range forces 

Massless scalars can also be stabilized by imposing 
both supersymmetry and chiral symmetry - the chiral 
symmetry enforces the masslessness of a fermion species, 
and the supersymmetry requires the scalar and fermion 
to be degenerate. Massless spin-two bosons (gravitons) 
can be stabilized by general covariance. If we could find 
spin models in which supersymmetry and/or general co- 
variance are emergent and robust with respect to generic 
local perturbations of the microscopic Hamiltonian, that 
would be very interesting from the perspective of funda¬ 
mental physics, apart from the potential applications to 
self-correcting quantum memory! 

Putting aside the microscopic origin and perturbative 


A. Quasi-topological phases 

In mathematical models of topologically ordered sys¬ 
tems, we usually assume the system has an energy 
gap separating the topologically degenerate ground state 
from the rest of the energy spectrum. Then we can imag¬ 
ine integrating out all quasiparticle excitations to obtain 
an effective theory with no propagating excitations, a 
topological quantum field theory (TQFT) which accu¬ 
rately describes the physics of the system up to correc¬ 
tions which are exponentially small in the system size. 
But strictly speaking this picture does not necessarily ap¬ 
ply to the topologically ordered systems observed in the 
laboratory, which may support gapless photon or phonon 
excitations. Furthermore, disordered systems may have 
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many low-lying localized bulk quasiparticle excitations 
with energy well below the topological gap. 

Bonderson and Nayak proposed the term quasi- 
topological phase for a system that has some of the fea¬ 
tures of a topological phase, but also has gapless excita¬ 
tions and therefore cannot be precisely described at low 
energy by a TQFT [15]. For example, Laughlin states 
are quasi-topological; their fractionally charged quasi¬ 
particles couple to the massless photon, and correspond¬ 
ingly the topological degeneracy on the torus is split by 
an amount that scales polynomially rather than expo¬ 
nentially with the system size. Furthermore, the anyon 
braiding statistics is modified by nonuniversal dynami¬ 
cal phases which are not described by the TQFT. Bon¬ 
derson and Nayak also proposed the term strong quasi- 
topological phase for a system in which some sector of 
the theory decouples from the gapless excitations and 
can be accurately described by a TQFT. For example, 
the Moore-Read state has a strongly quasi-topological 
sector. Although the Ising anyons carry electric charges 
which produce non-topological exchange phases, the dis¬ 
tinguishable fusion channels for a system of many Ising 
anyons span a topologically protected Hilbert space well 
described by the Ising TQFT. 

The toric-boson model, the toric code coupled to an 
exactly massless scalar field as in Eq. Q, is a strong 
quasi-topological phase. One might have thought that, 
just as distantly separated anyons have a potential en¬ 
ergy scaling algebraically with their separation, so large 
punctures on the lattice, which we may think of as fat¬ 
tened anyons, should also have an algebraically decaying 
interaction, and corresponding the splitting of the de¬ 
generacy on the torus should be algebraically decaying 
as for Laughlin states. But in the exactly solvable toric- 
boson model, the topological degeneracy on the torus is 
exact. Deforming the model with a small local pertur¬ 
bation might conceivably destroy the topological order 
as discussed in Sec. |Vj but in the “weak magnetic field” 
regime of Eq. (26) and (27) the perturbation may pro¬ 
duce only an exponentially small splitting of the ground 
state degeneracy, just as for the toric code without long- 
range interactions. 


B. Topological order at low temperature in a Unite 
system 

Our focus throughout this paper has been on the scal¬ 
ing of the memory time with system size, and in particu¬ 
lar on whether arbitrarily long memory times are achiev¬ 
able as a matter of principle. But we should emphasize 
that a two-dimensional topologically ordered system of 
fixed size could be a useful quantum memory if the tem¬ 
perature T is far below the topological gap A. 

The signatures of topological order, such as quantized 
Hall conductance and nontrivial topological entangle¬ 
ment entropy, disappear on an L x L lattice for system 


size larger than L 0 (T), where JTT, 42 

logLoCT)« A/T. (31) 

For L > L 0 (T) the system is likely to contain propa¬ 
gating thermally excited anyons, causing rapid decay of 
stored quantum information. But just as the quantized 
Hall conductance is a robust physical phenomenon, even 
though it disappears in the thermodynamic limit at any 
nonzero temperature, so a topological quantum memory 
with storage time enhanced by the inverse Boltzmann 
factor e A / T rnay be a valuable resource, even though it 
fails to meet our criteria for self correction. 


C. Self correction and fast propagation 

The Lieb-Robinson bound, which applies to any local 
Hamiltonian which is a sum of finite-norm terms, asserts 
that information propagates at finite speed. This bound 
is violated by the toric-boson model Eq. ([3|; therefore 
this model cannot arise from an underlying model whose 
Hamiltonian is a sum of terms, each with bounded norm 
and bounded range. 

When the toric-boson model on a L x L lattice is per¬ 
turbed by a suitable local Hamiltonian, anyons can prop¬ 
agate at speed 0{L). Consider dividing the lattice into 
two domains separated by a horizontal domain wall, as 
in Fig. [I] The lower half A is in the phase where plaque- 
ttes take the value W r = +1, and the upper half B is in 
the phase with W r = —1. At time t = 0, two anyonic 
excitations (each with W r = +1) are created deep within 
the B phase, and at the same time magnetic fields are 
applied to spins on a vertical line which extends from the 
lower anyon to the domain wall, and beyond into the A 
phase. This anyon is subjected to a long-range force at¬ 
tracting it toward the A phase. Were no magnetic field 
applied, the configuration with an anyon pair would be 
an exact energy eigenstate and would not evolve. But 
the applied magnetic field gives the lower anyon a finite 
effective mass, so it can respond to the attractive force 
by falling vertically. By the time it reaches the domain 
wall, this falling anyon has accelerated to speed O(L). 
Thus anyons can travel at unbounded speed in the limit 
of infinite system size, showing that the Lieb-Robinson 
bound does not apply. When the anyon reaches the do¬ 
main wall, it plunges through, becoming an anyon with 
W r = — 1 in the A phase. Forces attracting it back to¬ 
ward the B phase slow the anyon down until it comes to 
rest deep with the A phase, and then proceeds to oscillate 
between positions in the A and B phases. 

Since the divergent anyon chemical potential of the 
toric-boson model is directly related to the long-range 
forces which can accelerate anyons to unbounded speed, 
one wonders whether this connection holds more gener¬ 
ally. Generalized Lieb-Robinson bounds have been de¬ 
rived for systems which have long-range interactions scal¬ 
ing with separation r as r -a , showing that in D dimen- 
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excitations 



FIG. 1: A thought experiment to demonstrate the violation 
of the Lieb-Robinson bound in the toric-boson model. The A 
phase with W r = +1 is separated by a horizontal domain wall 
from the B phase with W r — — 1. An anyon pair is created 
at time t = 0 deep within the B phase, and at the same 
time a magnetic field turns on along a vertical line extending 
from one anyon to the A phase. This anyon, attracted by the 
A phase, falls vertically, reaching speed O(L) by the time it 
reaches the domain wall. 

sions the time required for information to propagate a 
distance r is bounded below by a logarithm of r for 
a > D [SB] and by a sublinear power of r for a > 2D [ST]. 
No limitations on propagation speed have been derived 
for a < D, the regime in which long-range interactions 
can cause the chemical potential to diverge. Since the in¬ 
finite anyon chemical potential results from the response 
of the topological medium arbitrarily far away, it’s tempt¬ 
ing to assert that a divergent chemical potential can oc¬ 
cur only in a system with constant propagation time in¬ 
dependent of r, though we do not have a rigorous general 
argument for that conclusion. 

Not only bosons, but also fermions can induce long- 
range interactions when fermion wave functions are de¬ 
localized. A familiar example is the RKKY interaction 
between nuclear magnetic dipole moments induced by 
exchange of conduction electrons, leading to an effective 
interaction strength proportional to 

~ —r[2kr cos(2fcr) — sin(2fcr)l (32) 

r 4 

where r is the separation of the dipoles and k is the con¬ 
duction electron’s wave number. Even though such effec¬ 
tive long-range forces can arise from fermion exchange, 
these forces cannot produce an infinite propagation speed 
or an infinite anyon chemical potential. The microscopic 
Hamiltonian is a sum of terms with bounded range and 
bounded norm, so the Lieb-Robinson bound implies a 
finite propagation speed. Furthermore, the energy cost 
of creating an isolated anyon is finite, because only a 
constant number of terms in the Hamiltonian, each with 
constant norm, are sensitive to the presence of the exci¬ 
tation. Indeed, this remark applies not just to fermions 
but also to spin models where the Hamiltonian is a sum 


of terms with bounded range and bounded norm. 


D. Self correction as a phase transition 

We expect the stored quantum information in a quan¬ 
tum memory to decay quickly at sufficiently high tem¬ 
perature. Therefore, a system which is self correcting 
is expected to undergo a phase transition at a critical 
temperature. In the low-temperature phase the mem¬ 
ory time increases without bound as the system size in¬ 
creases, and in the high-temperature phase the memory 
time is bounded above by a constant independent of sys¬ 
tem size. If the system has no phase transition at nonzero 
temperature, then we do not expect it to be self correct¬ 
ing at low temperature. 

For example, though Haah’s cubic code [20] and Mich- 
nicki’s welded code :22] [23] both have divergent energy 
barriers, in both cases the partition function can be com¬ 
puted exactly, and one finds that singularities occur only 
at zero temperature. The absence of a phase transition at 
nonzero temperature is compatible with the observation 
that the memory time is constant at nonzero tempera¬ 
ture, despite the increasing energy barrier. In contrast, 
the four-dimensional toric code, which is self correcting, 
has a phase transition at nonzero temperature associated 
with condensation of string excitations [52] . 

In fact the connection between self correction and 
phase transitions can be rather subtle, as Hastings et al. 
observed [53]. They noticed that in the toric code in six 
or more dimensions, the critical temperature T c associ¬ 
ated with thermodynamic singularities is higher than the 
percolation temperature T p , above which thermally ex¬ 
cited extended defects can have infinite size with nonneg- 
ligible probability; furthermore, the separation between 
the two temperatures grows parametrically with increas¬ 
ing dimension D. A related phenomenon occurs if we 
consider the toric code in D = 4 dimensions, but where 
the variables residing on lattice plaquettes take values in 
Zn rather than Z 2 . As has been known for some time, 
for N > 5 the four-dinrensional Zjy gauge theory has 
two phase transitions as the temperature varies [52j . In 
the low-temperature phase, defects are typically small 
and dilute, and the memory time increases exponentially 
with system size. In the high-temperature phase, stored 
information decays in constant time. 

To gain insight into the memory time in the 
intermediate-temperature region, it is helpful to consider 
classical spin systems storing classical information which 
also exhibit defect percolation below the critical temper¬ 
ature. For example, for the JV-state Potts model [53] 
with N > 5, the defect percolation temperature T p , 
above which thermally fluctuating domain walls have un¬ 
bounded size, is less than the critical temperature T c , 
above which the system does not spontaneously magne¬ 
tize. In the intermediate-temperature phase, although 
defects percolate, the free-energy barrier for logical er¬ 
rors is O(logL), and correspondingly the memory time 
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scales polynomially with system size [5'5j . This physical 
picture suggests that the memory time may also scale 
polynomially for high-dimensional toric codes, when the 
temperature is in between the percolation temperature 
and the critical temperature. Thus, while T < T c may 
be necessary for an exponentially scaling memory time, 
it does not seem to be sufficient in general. 

Hastings et al. also suggested that high-dimensional 
toric codes may be self correcting even in the superheated 
regime above the critical temperature, due to hysteretic 
behavior [53]. Mean-held theory predicts that the Potts 
model is a good classical memory even for T > T c [55] , 
but this mean-held prediction is not supported by nu¬ 
merical simulations. The memory time of the (slightly) 
superheated two-dimensional Ising model increases poly¬ 
nomially with system size L up to an optimal size L* 
comparable to the correlation length, but the IV-state 
Potts model has a first-order thermal phase transition 
for N > 5, and its stored classical information decays 
rapidly when the temperature is just slightly above the 
critical temperature. 

Another connection between self correction and phase 
transitions can be obtained from the analysis of master 
equations describing thermally fluctuating systems. For 
any classical spin system subject to Glauber dynamics, in 
which all correlations functions decay exponentially with 
distance, if the Gibbs state is the unique hxed point then 
the convergence to the Gibbs state is rapid m- Here 
“rapid” means scaling polynomially with system size; for 
some classical spin systems stronger results can be de¬ 
rived |58| , and in any case the classical memory time may 
in some cases be much shorter than the convergence time 
to the Gibbs state. Since exponentially decaying corre¬ 
lations are a typical feature of high-temperature phases, 
we may expect rapid mixing to the Gibbs state above the 
critical temperature, and therefore rapid decay of stored 
classical information. 


E. Perturbative origin of long-range forces 


In Sec. IV we discussed how long-range interactions 
between plaquette operators might enhance the quantum 
memory time of a two-dimensional topologically ordered 
system. But until now we have not considered in de¬ 
tail how such long-range forces might arise from a micro¬ 
scopic Hamiltonian which is a sum of terms with bounded 
range and bounded norm. 

Before we discuss the microscopic origin of long-range 
interactions, let’s first recall how topological order can 
arise in a two-dimensional spin system. The toric 
code provides a beautiful and instructive example of a 
frustration-free commuting Hamiltonian with a topologi¬ 
cally ordered ground state. But each term in the Hamil¬ 
tonian acts nontrivially on four qubits, while the nat¬ 
urally occurring interactions in a spin Hamiltonian are 
typically two-body terms. In fact, frustration-free com¬ 
muting models realizing any doubled (non-chiral) anyon 


model can be constructed [5H], but only by including 
even higher-weight many-body terms in the Hamiltonian. 
What spin models with physically plausible two-local in¬ 
teractions will exhibit topological order? 

It is known that the ground state of any two-body 
frustration-free commuting Hamiltonian has only short- 
range entanglement, and therefore cannot be topolog¬ 
ically ordered [SOI El]. This conclusion also applies 
to frustration-free commuting Hamiltonians with three- 
local interactions among qubits on any graph, or for 
three-local interactions among qutrits on a nearly Eu¬ 
clidean lattice [B2j. In this sense the toric code, with 
four-local interactions among qubits on a square lattice, 
is the optimal model of two-dimensional topological order 
based on an exactly solvable commuting Hamiltonian. 

Therefore, to obtain a topologically ordered ground 
state from a two-local Hamiltonian, we must be willing 
to consider noncommuting Hamiltonians, which may not 
be exactly solvable. One approach to constructing such 
models is the method of perturbative gadgets [63], 641 . 
The idea is to approximate a target Hamiltonian Ho, 
which is commuting and has a topologically ordered 
ground state, with a two-local Hamiltonian H. We may 
express H as 


H = H 0 + V; (33) 

because the topological order of H 0 is perturbatively sta¬ 
ble, the quantum memory properties of the two-local 
Hamiltonian H will mimic those of the target Hamil¬ 
tonian H 0 if the perturbation V is a sum of sufficiently 
small local terms. 

But can we also obtain long-range plaquette-plaquette 
interactions of the form — ]T] r r , J r yW r W r ' as a pertur¬ 
bative approximation to a two-body Hamiltonian? On 
a square lattice, each term in this Hamiltonian involves 
eight lattice sites, and the number of terms scales like the 
square of the system’s area, where each term involves the 
exchange of at least one gadget particle, so the “density” 
of gadget particles per unit area is divergent. We have 
not found a way for the gadget particles to couple to the 
plaquette variables in the desired way, or for the nonlo¬ 
cal Hamiltonian to be perturbatively close to a two-local 
one. 

In fact, as already noted in Sec. |VIIC| we do not expect 
the infinite anyon chemical potential of the toric-boson 
model to be realizable by any Hamiltonian which is a sum 
of terms with bounded range and bounded norm, simply 
because the energy cost of a localized excitation should 
be finite if the excitation occupies a finite lattice volume. 
Or at any rate, the toric-boson model could provide a 
good description only in a scaling limit where the lattice 
spacing becomes small compared to the anyon size. 

A further potential problem is that long-range interac¬ 
tions, whatever their origin, might interfere with the per¬ 
turbative stability of the target Hamiltonian, and hence 
invalidate the perturbative gadget method for achieving 
topological order. Arguments for the perturbative sta- 
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bility of topological order formulated in [0G2I do not 
apply if the perturbation V includes terms which decay 
algebraically (rather than exponentially) with distance. 
Indeed, systems with long-range interactions can often be 
well described using mean-field approximations in which 
entanglement is assumed to be short range. If a topologi¬ 
cally ordered system with long-range plaquette-plaquette 
interactions can arise from a two-local Hamiltonian at all, 
approximation schemes going beyond conventional per¬ 
turbation theory may be needed to analyze the stability 
of the system’s stored quantum information [55] . 
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